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The lattice provides a powerful tool to non-perturbatively investigate strongly coupled super- 
symmetric Yang-Mills (SYM) theories. The pure SU(2) SYM theory with one supercharge is 
simulated on large lattices with small Majorana gluino masses down to about antg = 0.068 with 
lattice spacing a ~ 0.125 fm. The gluino dynamics is simulated by the Two-Step Multi-Boson 
(TSMB) and the Two-Step Polynomial Hybrid Monte Carlo (TS-PHMC) algorithms. Supersym- 
metry (SUSY) is broken explicitly by the lattice and the Wilson term and softly by the presence of 
a non-vanishing gluino mass. However, the recovery of SUSY is expected in the infinite volume 
continuum limit by tuning the bare parameters to the SUSY point in the parameter space. This 
scenario is studied by the determination of the low-energy mass spectrum and by means of lattice 
SUSY Ward-Identities (Wis). 
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1. Introduction 

Non-perturbative dynamical effects like confinement and the spontaneous breaking of a dis- 
crete chiral symmetry are expected to occur in the pure JV = 1 supersymmetric Yang-Mills (SYM) 
theory with SU(A'^c) as the color gauge group. The confinement is realized by colorless bound states, 
while the anomaly leaves only a discrete subgroup of the chiral transformations invariant, which is 
spontaneously broken by the non-zero gluino condensate. Continuum low-energy effective actions 
[§, ^] predict two chiral supermultiplets at the lower-end of the particle spectrum, each containing 
two scalars with opposite parities and one Majorana fermion. A massive gluino introduces a soft 
breaking of supersymmetry (SUSY) resolving the degeneracy in the supermultiplets. 

We continue here the investigations of the DESY-Munster collaboration on non-perturbative 
features of the = 1 SU(2) SYM theory using lattice simulations (for reviews, see [|] and ref- 
erences therein.) The Wilson formulation for gluino action is applied while for the gauge (gluon) 
action we used pure and improved Wilson lattice actions. The Wilson formulation for the gluino 
action breaks both SUSY and chiral symmetry explicitly, and the presence of a non-vanishing mass 
term breaks SUSY softly, too. The recovery of SUSY and chiral symmetry is expected in the in- 
finite volume continuum limit by an appropriate fine tuning of the bare gauge coupling and the 
bare gluino mass to their critical values. To this end, the mass spectrum of low-lying bound states 
of the theory is computed. Furthermore, the lattice SUSY Ward-Identities (Wis) can be checked 
independently, indicating that SUSY is restored in the continuum limit. 

2. Simulation details 

Table 1: Algorithmic parameters for the j3 — 2.3 , 16^ -32 TSMB runs. A^vc the acceptance rate in the 
noisy correction step, x^^"^ is the autocorrelation times of the plaquette. For e, X, ni, «2 see text. 



Run 
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# Sweep 


Anc % 
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ni 


«2 


(«) 


0.1955 


12500 


50-80 


167.6 


2.0 


10"^ 


4.0 


40 


800 


ib) 


0.1960 


23500 


50-80 


181.1 


4.0 


10"^ 


4.0 


40 


1800 


(c) 


0.1965 


18000 


50-62 


254.2 


4.0 


10-^ 


4.0 


40 


1800 



To simulate gluinos on the lattice, we used the Curci-Veneziano (CV) effective gauge action 
[^. A range of lattices were simulated using the Wilson formulation for both the gluino and the 
gluon parts of the lattice action. In recent simulations, stout-smearing of the gauge fields used in 
the Dirac-Wilson operator has been applied to reduce the fluctuations of the smallest eigenvalues 
of the fermion matrix resulting in only a few exceptional configurations. Furthermore, in the gauge 
sector we used an improved lattice action, namely the tree level improved Symanzik (tlSym) gauge 
action (For more details see ref. [Q]). 

We used two variants of dynamical updating algorithms: 1) at gauge coupling j8 = 2.3 three 
ensembles at hopping parameter values K = 0.1955, K = 0.196 and K = 0.1965 were prepared by 
the Two-Step Multi-Boson (TSMB) algorithm |§] on 16^^ • 32 lattices. 2) at gauge couphng jS = 2.1 
several 16^ • 32 and 24^ -48 lattices were generated by the Two-Step Polynomial Hybrid Monte 
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SYM on 24c48, TISym 5=1.60 k = 0.1 570 



<P>=0.625970(24), :,n,=7.2 




Figure 1: History of the average plaquette and the smallest eigenvalues. 



Carlo (TS-PHMC) algorithm The input parameters of the TSMB runs are summarized in Table 
[T] and the TS-PHMC runs input parameters can be found in Table 1 of Ref . [Q] . These algorithms are 
based on the multi-boson representation of the fermion determinant by approximating the inverse 
fermion matrix by a polynomial of order n\ in an interval [£, A] and a stochastic noisy correction 
(NC) step realized by a second polynomial of higher order ?i2- The history of the average plaquette 
and the smallest eigenvalue of one run using stout-smearing is displayed in Fig. |l|. 



3. Confinement and physical scale 

Analogy with QCD suggests that the theory presents confinement with color singlet particle 
states. Confinement can be studied by computing the static quark potential Vif) as a function of 
the spatial separation r between sources in the fundamental representation of the gauge group. The 
basic quantities to determine the static potential are the Wilson loops. The static potential Vir) in 
Fig. m was extracted from the single-exponential fit as described in [|l|]. We applied APE smearing 
(£ap£ = 0.5, Nape = 4) for Wilson loops to improve the signal and to reduce the contamination 
with excited states. 

The physical scale is set by the value of the string tension or the Sommer scale parameter ro/a. 
These are determined by two independent methods: using the force method based on the Creutz- 
ratios and by fitting V (r) directly. At lighter gluino mass we obtained ro/a ~ 8 for TSMB runs 
and ro/<3 ~ 4 for the TS-PHMC runs. Using = 0.5 fm, the lattice spacings are a ^ 0.06 fm and 
a « 0.125 fm for the TSMB and TS-PHMC runs, respectively. The physical size of the simulated 
boxes is then ~ (Ifm)^ for the TSMB sample and ~ (2 - Sfm)^ for the TS-PHMC sample. 



4. Supersymmetry and chiral limit 

The gluino mass can be determined by means of the lattice SUSY Ward-Identities (Wis). Since 
it is assumed to behave linearly in 1/k, one can extrapolate the inverse hopping parameter to the 
massless gluino limit, where both SUSY and chiral symmetry emerge. In that way the critical value 
of the hopping parameter can be determined. However, a second independent method based on 
OZI arguments [Q] can be used to estimate K^r, where the adjoint pion a — n becomes massless 
in the chiral limit. We observed good agreement between Kct estimated from these two methods. 
In the case of the TSMB ensembles we found ~ 0.1969, and for the TS-PHMC ensembles 
Kcr — 0.2033. The linear extrapolations for both methods are shown in Fig. ^ and Fig. ^ 
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Lattice: 16'. 32 P=1.6 K=0.2 PHMC-tlsym 




Figure 2: The static potential. 




5. Spectrum results and conclusion 

In the SUSY point the low-energy effective actions describe the lower-end particle spectrum 
by two massive chiral supermultiplets ||]. The predicted composite states are color singlet bound 
states: two scalars (a — fo, 0+ glueball), two pseudoscalars (a — rj', glueball) and two spin-^ 
states X (gluino-glueballs). The physical states of the two supermultiplets result in general from 
mixing of these states. In lattice simulations correlation functions of interpolating fields with the 
corresponding quantum numbers are computed for large time separations such that the lightest state 
with these quantum numbers (unmixed state) dominates. By this method, the masses of the bound 
states are extracted. The results for the mass spectrum are presented in Fig. ^ as a function of the 
inverse hopping parameter. The vertical line indicates the massless gluino region. 

The results of TSMB simulations in a (Ifm)^ physical volume indicate the degeneracy of 
a — rj' and x masses and hence the two bound states may occur in a supermultiplet. However, in the 
larger volumes of TS-PHMC runs a mass difference between these two states has been observed, 
which appears to be larger than the expected value of the mass-splitting from the soft breaking, 
given by the gluino mass. We conclude that the TSMB results are dominated by large finite size 
effects. Furthermore, the masses of three bound states a — /o and 0+ glueball converge to one 
point near the region where SUSY is expected to emerge, and hence this is a signal of the heavier 
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5,1 



Figure 4: The gluino mass (left) and the pion mass square (right) as function of the inverse hopping param- 
eter (TS-PHMC 16^^ ■ 32). 
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Figure 5: Low-lying spectrum of JV = 1 SU(2) SYM theory. The mass in lattice unit as function of the 
bare gluino mass 1 / fc. 



supermultiplet, while the a — rj' remains lighter and it could belong to the lightest supermultiplet. 
Simulations on a finer lattice will clarify the impact on these findings of explicit SUSY breaking. 

Our computations were performed on the Blue Gene L/P and JuMP at the Neumann Institute 
for Computing (NIC), the Cluster of the University of Miinster and the Cluster of RWTH Aachen. 
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